SVZ © l/q 2 expansion versus some QCD holographic Models 



F. Jugeau 3 , S. Narison b ' 1,2 , H. Ratsimbarison c 



"Instituto de Fisica, Universidade Federal do Rio de Janeiro, Caixa Postal 68528, RJ 21941-972, Rio de Janeiro, Brazil. 
b Laboratoire Particules et Univers de Montpellier, CNRS-IN2P3, Case 070, Place Eugene Bataillon, 34095 - Montpellier, France. 
1 Institute of High-Energy Physics of Madagascar (iHEP-MAD), University of Antananarivo, Madagascar 



Abstract 

. . . ( Considering the classical two-point correlators built from (axial)-vector, scalar qq and gluonium currents, we confront results 
obtained using the SVZ © l/q 2 expansion to the ones from some QCD holographic models in the Euclidian region and with 
_| negative dilaton 4>,(z) = -|c 2 |z 2 . We conclude that the presence of the l/^-term in the SVZ-expansion due to a tachyonic gluon 
mass appears naturally in the Minimum Soft Wall (MSW) and the Gauge/String Dual (GSD) models which can also reproduce 
^ ]semi-quantitatively some of the higher dimension condensate contributions appearing in the OPE. The Hard- Wall model shows 
C$ a large departure from the SVZ © l/q 2 expansion in the vector, scalar and gluonium channels due to the absence of any power 
^) corrections. The equivalence of the MSW and GSD models is manifest in the vector channel through the relation of the dilaton 
parameter with the tachyonic gluon mass. For approximately reproducing the phenomenological values of the dimension d — 4, 6 
OS condensates, the holographic models require a tachyonic gluon mass (a s /n)A 2 m -(0.12 ~ 0.14) GeV 2 , which is about twice the 

t ! fitted phenomenological value from e + e~ data. The relation of the inverse length parameter c, to the tachyonic gluon mass also 

<~| shows that c, is channel dependent but not universal for a given holographic model. Using the MSW model and M p = 0.78 GeV as 
Q^idnput, we predict a scalar qq mass Ms ~ (0.95 ~ 1.10) GeV and a scalar gluonium mass Mq » (1.1 ~ 1.3) GeV. 

^j^'Keywords: QCD spectral sum rules, AdS/QCD, Conformal Field Theory. 
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1. Introduction 

The SVZ-sum rules approach [1, 2] have been used suc- 
cessfully for understanding the properties of the lowest mass 
hadrons in QCD. It is based on a duality between the measured 
spectral function and the semi-perturbative QCD expression us- 
ing the Operator Product Expansion (OPE) in terms of the quark 
and gluon condensates or their mixing. These condensates are 
assumed to give a good description of confinement at moder- 
ate energies. More recently [3, 4], it has been also shown that 
UV renormalon contributions can be mimics by a 1 /q ,2 -term in- 
duced by a tachyonic gluon mass while its duality with these 
large order terms of the PT series has been studied in [5]. 
On the other hand, QCD holographic models have been re- 
cently extensively studied following the ideas in [6, 7] that one 
can have a duality between the large N c limit of a maximally 
N-A superconformal SU(N C ) gauge theory in n dimensions 
and the supergravity limit of a superstring/ M-theory living on 
a (n + 1) Anti de Sitter (AdS) space ® a compact manifold. 
Such a gauge/gravity duality has been applied to QCD by map- 
ping it onto gravitational effective theories in 5-dimensions on 
an AdS background. Different forms of the gauge/string corre- 
spondence for studying the hadron spectra in Minkowski space 
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have been proposed in the literature. In the so-called top-down 
approximation, one tries to keep the underlying string structure. 
In so doing, one usually starts with a theory in 10-dimensions. 
Its reduction to 5-dimensions in the supergravity limit contains 
in general additional higher derivatives which are 1 / y[N^ (or 
the Regge slope a') corrections. A model of this approach has 
been discussed in [8] and used in different phenomenological 
applications in [9-1 1]. 

In this paper, we shall consider the so-called Soft- and Hard- 
Wall models [12, 13] of the bottom-up approach and the sim- 
plified version of the Gauge/String Duality (GSD) [8] used in 
[9-1 1] . We compare the results obtained for the two-point func- 
tions from these holographic models with the SVZ © l/q 2 ex- 
pansion. 

In so doing, we consider the two-point correlators: 



U(q 2 ) 



-J 



d A xe lqx {Q\T{J{x), /(0)}|0>, 



built from the local (axial-) vector current: 

or the local scalar quark current: 

J A (x) = 4>T A iKx), 

where i/s is the light quark field and T A = (T° — 1 / J2n /, T") 
with T" = A a /2 the S U(3) f generators (a = 1, . . . , 8) and the 
normalization tr(T A T B ) = 6 AB /2. We shall also consider the 



(1) 



(2) 



(3) 



gluonium correlator built from gluonium field strength: 
J G (x) = /3(a s )G%G» va (x) : /3(a s ) = # (^) + • ■ 



(4) 
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with B\ = -(1/2)(11 - 2n//3) for nj active flavours. We 
shall evaluate the previous two-point correlators using either 
the previous holographic models or the SVZ expansion [1, 2] 
including quark and gluon condensates and the l/q 2 contribu- 
tion from tachyonic gluon mass which goes beyond the stan- 
dard SVZ-expansion [3]. The latter term has been motivated 
to phenomenologically modelize the UV-renormalon contribu- 
tions and has been shown in [5] to be dual to large order PT 
series. 

2. The vector two-point function II {q ) 

The two-point correlator has the Lorentz decomposition: 
Tl^ ah (q 2 ) = -8 ab (g" v q 2 - q"q v ) IT^ + S^^Ilf (5) 

where the upper indices 1 and refer to the spin 1 (transverse) 
and spin (longitudinal) hadronic components. 

( 1) 9 9 

• Tly (q ) using the SVZ expansion © l/q -term 

One can evaluate the two-point correlation function using the 
SVZ expansion ffi l/^-term in the Euclidian region. One ob- 
tains to leading order in a s and using the normalization of the 
current in Eq. (2) [1-3]: 

(^)n<V> - + 

ln\{a,G 2 ) 896 3 pa,(i//i//} 2 
+ \3)- 



q-* 81 q u 

-q 2 3607 \ (a s G 2 ) 2 
-I 531 ° g — + l22-] 259V' 



(6) 



bulk action for a gauge vector field over the bulk-volume *V is 
given by [13]: 



- f d^Ve-^—trF 2 



(10) 



where k is the 5 -dimensional Newton constant and F" V ^ MN = 
3mV^ - d^V a M is the abelian part of the field strength tensor. 
The corresponding equation of motion reads: 



dui^ge-^Ff^) 







(11) 



which, in terms of the Ad Fourier-transformed field V a M (q, z) of 
y a M {x, z), becomes in the axial-like gauge V = 0: 

' R _ rT ,f,i „ - . \ 2 ^ 



d z ^-e-^ ) d z V^q,z)yq 2 -e- <m V' la (q,z) = (12) 



and: 



q lt V tia (q,z) = 0. 



(13) 



The equation of motion of the bulk-to-boundary propagator 
reads: 



2clz + 



l)d z -q 2 



V\%,c 2 ,z 2 







(14) 



with the general solution (z > 0): 



2 2 

V(-c vZ 



B(q 2 )c 2 ,z 2 iF l [^ + l;2;c 2 ,z 2 \( 



where 



a = r i + 



4c 2 



(16) 



where we have used the contribution of the d — 8 condensates 
within the vacuum saturation assumption but including the 1 /N c 
corrections from [14]. A 2 is the tachyonic gluon mass with a 
phenomenological value [3, 15]: 



i^A 2 = -(0.06 ~ 0.07) GeV 2 , 
and with an absolute approximate upper limit: 
(^■)|/1 2 | < (0.12 ~ 0.14) GeV 2 . 



(7) 



(8) 
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• Tly (q < ) from the MSW model with a positive dilaton field 

The vector two-point function in the Minimal Soft Wall (MSW) 
model has been studied in [13, 16], where the case of a 
Minkowskian metric and a positive background dilaton field 
<l) v ,(z) = c 2 z 2 with c 2 > has been considered. The models 
are characterized by the metric: 



ds 2 = e 2A(z} [rj liv dxf l dx v + dz 2 ) = g M N<ix M dx N 



(9) 



where A(z) = - ln(z/R) for pure AdSs with radius R; rj^ v = g MV 
is the 4D Minkowski metric tensor with signature (-+++). The 



and B are integration constants with respect to z; U is the Tri- 
comi confluent hypergeometric function and i F\ is the Kummer 
confluent hypergeometric function. In the UV limit z — > 0, this 
solution reads: 



i + 



2 2 
q Z 



4 

+2y E - 1 



i(c 2 z 2 ) + ( a(£t + 1 



+ B{q 2 )c 2 v z 2 + 0{z") , (17) 



where the digamma function iff(a) has an infinite set of simple 
poles located at a = —n (n 6 N) with residues -1 and 



1 V" 1 Blk 1 

tl/(a) -\na > — 

* 2a f-f 2k a 2k 



(18) 



where B2 = 1/6, B\ = -1/30,... are Bernoulli numbers. In the 
Fourier g-space, the two-point correlator is obtained by means 
of twice the functional derivative of the critical bulk action with 
respect to the boundary field Vq: 



-^(/Tr-^V'KV) 



-i(2n) 4 f 



d 4 k^ 



V a =0 ' 



(19) 



from which one can deduce: 



4V) 



io 20) 



• Yly {q ) from the MSW model with a negative dilaton field 

One can do an analysis similar to the case of the positive dilaton 
field [17]. Defining the negative dilaton field by: 



where the critical field V cr j t (q,z), also solution of Eq. (14), is 
related to its restriction Vo{q) on the boundary space by means 
of the bulk-to-boundary propagator V(K , c 2 z 2 ): 



(21) 



Using the solution in Eq.(20), one can deduce in the Minkowski 
space [13]: 



(kg 



2\ 



R 



n ( ' ) = --log^-^ + ^-^^ + 0( 3 -|(22) 



q 



2 ct 16 cl 



q 2 3 q 4 15 q* 



1 



14 



where we have dropped constant terms which are physically 
irrelevant; v is the usual subtraction constant. It is impor- 
tant to notice that due to an "accident" numerical cancellation 
among different terms, the coefficients of c v , c l v ° and c 12 van- 
ish in the expansion. By comparing this result obtained in the 
Minkowskian region with the SVZ one in the Euclidian space, 
one should first Wick rotate the action in Eq. (10) such that 
q 2 M > becomes -q\ > (in the signature (-++++) and where 
the subscripts M and E denote respectively Minkowski and Eu- 
clidian 3 ) and change q 2 M to —q 2 E '. 



(23) 
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Tly (q ) from an extension of the MSW model 



The absence of the l/q 6 term in the MSW model due to some 
miraculous cancellation of its coefficients shows its departure 
from the SVZ-expansion. One can try to improve the model by 
adding the contribution of the large q 2 limit B of the B-term in 
Eq. (15). In this way, one obtains the additionnal contribution: 



R 



6U v (q 2 ) = ~B(q 2 )c 2 v 
1 



(24) 



where B(q 2 ) is the asymptotic form of B(q 2 ) in the large q 2 
limit. The presence of this term is however unlikely since this 
solution gives a divergent action S e J7 when z — > +°o contrary 
to the prescription in MSW model where the solutions corre- 
spond to a regular action. Keeping this solution in B would be 
similar to another choice of regularization prescription. In this 
way, one can obtain new terms: 



B(q 2 ) = B + B 2 ^ + B 4 ^ + B 6 ^- 

q^ q* q° 



(25) 



where B,- are arbitrary coefficients such that the modified model 
looses its predictivity and becomes less interesting. 



3 In the remaining part of the paper after Eq. (23), we shall omit these indices 
for simplifying the notation. 



<D v (z) = -\c 2 v \z 2 , 
the corresponding bulk propagator reads: 

y <0 = v(c?«-» |c v |V |CvlV 

= V{z -» 0) - |c v |V + 0(z 3 ) 



(26) 



(27) 



where V is the bulk propagator for a positive dilaton field de- 
fined in Eq. (15). Using Eq. (20) for V<o, one can deduce in the 
Minkowski space: 



n ( > 2 ) -ia) = n% 2 >c ?>0) + ^ 1 



(28) 



where there is an additional term (interpreted as a massless un- 
physical pole in the literature) compared with the expression 
> 0. Expanding the previous expression for 



of Ilf for c 2 
q 2 — > +oo, one obtains: 



R 



n 



id 



1 q 2 



q 2 + 3 q 4 



15 q s 



2\c 2 \ 



(29) 



where the additionnal term proper to c 2 < changes the sign 
of the total contribution of the l/q ,2 -term. Translating the re- 
sult to the Euclidian region by changing q 2 > to q 2 < 0, 
the final result becomes the same as in the case of a positive 
dilaton field in Eq. (22) when translated to the Euclidian re- 
gion. From this exercise, we show that the expression of ITV 
for <t>!,(z) = ±|c 2 |z 2 remains the same thanks to the existence of 
this additional 1 /g 2 -term, in the case of a negative c 2 , where its 
physical interpretation remains to be clarified. 

• n ( y \q 2 ) from Hard Wall (HW) model 

Evaluating the vector correlator in the HW model, one ob- 
tains [18]: 



-^|n ( ^V)--2log^-^ + ln2, 



1 



(30) 



where the perturbative part is the same as in the MSW model 
which is natural because the two models only differ in the IR 
region. However, it is rewarding that in the HW model, has 
no power corrections at all indicating its large departure from 
the SVZ-expansion. 

• Tly (q ) from a gauge/string dual ( GSD) model 

Using a gauge/string dual (GSD) model [8], n^ 1 ' has been eval- 
uated until the 1 /q 2 term in [9] in the background given by: 

R 2 

ds 2 = —e--- CAZ ~ (j^tfcW + dz 2 ) = g MN dx u dx N , (31) 



with the result: 



n ( v V) 



log ^ 



3ca 
q 2 



(32) 



The coefficient of the l/q 2 term can be compared with the one 
from the MSW model in Eq. (22) from which one can deduce: 



2c v = 3c A 



(33) 



By using a Weyl tranformation, one can show that the MSW 
model and GSD model are indeed equivalent in the vector chan- 
nel taking into account the constraint in Eq. (33) though [9] 
derives a slightly different constraint 2c 2 , — ca/2 which is due 
to the fact that he started from a string-like action in the back- 
ground (31) which differs from the one in Eq. (10). 

• SVZ © 1 /q 2 -expansion versus some QCD holographic models 

- Matching the SVZ © 1 /^-expansion result in the Euclidian 
region in Eq. (6), with the previous results from different holo- 
graphic models, one can remark that the HW model does not 
present any power corrections at all indicating its large depar- 
ture from the SVZ-idea. 

- From the previous analysis, one can notice that the quadratic 
correction introduced by [3] appears naturally in the MSW 
model and GSD model. This quadratic correction could be a 
remnant of stringy effect in QCD as suggested by Zakharov. 
Confronting the expressions in Eqs. (6) and (22) by doing the 
substitution in Eq. (23), leads to: 



kg] 



Yin 2 



and 2cf = 1.05 



(?)• 



(34) 



Using the previous value of A 2 in Eq. (7), one can deduce: 



c 2 , = -(0.03 ~ 0.04) GeV 2 < : 



(35) 



where the sign of c 2 is consistent with a tachyonic gluon. 
- Using this result, one can predict for the d — 4 condensate: 



(a s G 2 ) = -4 * (1 ~ 2) X 10~ 3 GeV 4 



which is relatively small compared to the phenomenological 
value [15, 19-26]: 



(a s G 2 ) = (0.07 + 0.02) GeV 4 , 



(37) 



and to the original SVZ estimate (a s G 2 ) * 0.04 GeV 4 [1] which 
has been recovered in [10] from a direct evaluation of the Wil- 
son loop in the GSD model. An analogous underestimate oc- 
curs for the d — 8 condensate obtained using a vacuum satu- 
ration assumption but including the l/N c corrections. In this 
d = 8 case, one can also notice that the MSW model does not 
reproduce the log —q 2 /v 2 -term appearing in the OPE. 
- Inversely, we can instead use the value of (a s G 2 ) in Eq. (37) 
to fix |c 2 |. We obtain in this way: 



\c 2 \ = 0.23 GeV 2 => j \A 2 \ = 0.45 GeV 2 



(38) 



Instead, one can also use the relation between |c 2 | and the p- 
meson spectrum [13, 17, 27]: 



M 2 _ =4|c 2 |(h + 1) 



(39) 



coming from the pole of the digamma function ifr (l + q 2 /4\c v \ 2 ^) 
appearing in Eq. (27) when q 2 = -M 2 , Demanding that the 
MSW model reproduces the lowest ground state p-meson (n = 
0), one can deduce: 



|c 2 U 0.15 GeV 2 



0.30 GeV . 



(40) 



The absolute values in Eqs. (38) and (40) are much larger than 
the phenomenological upper value of about (0.12 ~ 0.14) GeV 2 
obtained in [15]. Within the crude approximation where these 
numbers have been derived, we consider that the agreement 
with the phenomenological fit is remarkable. It is also impor- 
tant to notice that the matching of the MSW model results in the 
Euclidian region with the S VZ-expansion is essential for fixing 
the sign of c 2 which cannot be done from the alone fit of the 
spectrum in the Minkowski space due to the fact that the bulk 
potential from the Schrodinger-type equation behaves like c 4 . 

- As already mentioned earlier, the d = 6 condensate is absent 
in the MSW model due to some numerical cancellation of its 
coefficient. The extension of the MSW model can cure this dis- 
ease by the introduction of a non-zero integration constant B 
like in Eq. (25). However, within this framework, the MSW 
model looses its predictivity due to the arbitrariness of B and 
to the lack of convergence of the corresponding effective action 
S e ^J when z — > +°o. 

- A result similar to the one in the MSW model has been ob- 
tained by [9] in the GSD model. 



3. The IT, 1 ' (q 2 ) = H^Uq 2 ) - ET V) correlator 

LR V A 

YIlr controls the S U (h)lxS U(ii)r global chiral symmetry [2] 
and controls the assumption of the superconvergence of the 1 st 
(which involves the sum of the spin 1 and contributions) and 
the 2 nd (which involves the spin 1 or contribution) Weinberg 
sum rules for q 2 — » oo [28]. Solving these sum rules in the 
(36) chiral limit m q = and the current algebra relation in Ref. [29], 



Weinberg has derived the well-known relation: 
M Al - V2M p . 



(41) 



It has been shown within QCD [30] that the perturbative (PT) 
contribution due to the light quark masses breaks the 1 st Wein- 
berg sum rule to order a s m 2 while the 2 nd sum rule is broken to 
leading order m 2 log m 2 of PT. 



• n^(^ 2 ) in the SVZ® l/q 2 expansion 

In the chiral limit m q = 0, and using the SVZ expansion, one 
also finds a cancellation of contributions of terms of dimen- 
sion d < 4 due to m q {qq) and the flavour independent l/q 2 and 
(a s G 2 ) contributions. One obtains [1, 2, 31, 32]: 



q 6 Yl { ll = 87T P a s (uu) 2 

6tt(1+0) 



q°Yl 



32n 



-pa s {uu) 



(42) 



where p 2 indicates a deviation from vacuum saturation. An 
extraction of the four-quark condensate from e + e~ and r-decay 
data and the light baryon spectra gives [19, 20, 33]: 



pa s (uu) 2 (4.5 ± 0.3) X 10" 3 GeV 6 



Yl^iq 2 ) in the HW model 



One obtains in the Hard-Wall model [18]: 



An 2 



(43) 



where all contributions come from due to the fact that Yl v 
has no power corrections in this model. However, this feature 
is completely different from the SVZ-expansion. 



(44) 



in 



4. The scalar two-point function tt s (q 2 ) 

• Ylsiq 2 ) using the SVZ expansion © l/q 2 -term 
We pursue our investigation by evaluating the two-point func- 
tion of the scalar current using the SVZ © l/q 2 expansion. We 
shall work in the chiral limit m q = 0. Detailed derivations of 
the expression have been extensively discussed in the literature. 
To leading order in a s , the result is [1,2]: 



q 



log 



X 

«2 



n (a s G 2 ) 
~3~ 



log 



1 

„2 



llbn pa s {uu} 2 
~ ~27 q 1 
3677 \ {a s G 2 ) 2 
192 j 1728<? 6 ' 



(47) 



where we have used the result obtained by [14] for the d 
condensate contributions. 



Tl ( H(q 2 ) in the MSW model 



In the MSW model, the linear equation of motion of the axial- 
vector bulk fields (which describe holographic ally the axial res- 
onances) reads as: 

d z (-e-*®d z A x y qt-e-^A^ 

R -o (Z )^Vz) 2 



-A ± = . 



(45) 



v(z) is the function which describes the chiral symmetry break- 
ing. Typically, v(z) has two contributions: one which depends 
on the mass m q of the light quark while the other depends on the 
chiral condensate (qq). As a result, when m q = and keeping 
only the contributions of condensates d < 4, v{z) 2 vanishes and 
the vector and the axial sectors are completely equivalent [see 
Eq. (12)]. In the MSW model, the expression of Tim similar 
to the one in Eq. (44) has been derived in [34] from a relation 
of YIlr with the axial-vector-vector (AW) vertex following the 
evaluation of this vertex in [35]. We plan to come back to this 
analysis in a future publication. 



• SVZ © l/q -expansion versus holographic models 

One can compare the previous expressions (42) and (44) in the 
large N c limit. Using: 



(-) = — *- 

\nl -Pi 1 



:in c 



(46) 



one can notice that the previous expressions exhibit qualita- 
tively the same n 2 and large N c behaviours 4 but with different 
numerical values. 



• n s (q 2 ) from the MSW model 

The expression of this correlator is known in the literature for 
the case of a positive dilaton 5 . It reads in the Euclidean space 
[37, 38] 



I IW) = 



-(q 2 -2c 2 )log 



4 c° 28 c° 



56 c 



2c\_ 
3q 2 



(48) 



3 q 4 15 q 6 15 4 8 

where we have dropped constant terms which are physically 
irrelevant; v is the usual subtraction constant. We do the substi- 
tution in Eq. (23) for transfering this result into the Euclidian 
region, where one can also notice that the result remains the 
same for a negative dilaton. 

• SVZ © 1 j q 2 -expansion versus MSW model 

Matching the result from the two different approaches in the 

Euclidian region, one can deduce: 



R 



Nr. 



cl = 2\ 



(49) 



where c 2 is negative due to the tachyonic value of the gluon 
mass. Using this relation and the phenomenological value of 
a s A 2 in Eq. (7), one can predict: 



(a s G 2 ) = -c 4 s * (0.02 - 0.04) GeV 2 

7T 



(50) 



For the d = 6 contribution, the negative value of c 2 is crucial 
for matching the two approaches. One can deduce: 



paAqqf = * 1.1 X 10~ 3 , 



(51) 



compared with the phenomenological estimate in Eq. (43). 
One can notice that the AdS approach tends to underestimate 



4 The proper behaviour in the large- limit of the HW model has been e.g. 
discussed in [36]. 



5 One should notice that Ylsiq 2 ) does not have any short distance l/q 2 " 
power corrections in the HW model. 
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the value of the d — 4 and d — 6 condensates compared with 
the most recent phenomenological values. Alternatively, if one 
takes seriously the value of (a s G 2 ) = 0.07 GeV 4 from more 
recent phenomenology, one can predict: 



0.34 GeV 



(52) 



while a deviation by a factor 2 above the vacuum saturation of 
the four-quark condensates leads to: 



-0.30 GeV 



(53) 



where the sign has been fixed by the matching of the d — 6 
condensate contribution which requires again a tachyonic gluon 
mass. Therefore, our analysis shows that matching the MSW 
model approach with the S VZ expansion requires a value of the 
tachyonic gluon mass to be in the range: 



(*)■ 



-(0.15 ~ 0.17) GeV 



(54) 



where the absolute value is slightly higher but consistent with 
phenomenological determinations. Moreover, it is more im- 
portant to notice that the tachyonic origin of the gluon mass is 
necessary for matching the MSW model and the SVZ © l/q - 
expansion in the Euclidian region in order to have the correct 
signs of the d — 2 and d = 6 contributions. 

5. The scalar gluonium two-point function 

• Tlc(q 2 ) from SVZ® l/q 2 expansion 

Using SVZ-expansion for large q 2 in the Euclidian region and 
including the quadratic correction introduced by CNZ [3] for 
modelling UV renormalons and/or large order terms of the PT 
series, the expression of the two-point function reads for q 2 < 
and to leading order in a s [39, 40]: 



B(a s ) 



IW) 



(q 2 ) 2 



A 2 q 2 



log 



2n 2 



^ 3 G 3 > 



■{a s G 2 ) 



r 



2\2 



+ 9tt j (a s G z ) 



(55) 

8 a s 

where we have used a vacuum saturation including the 1 /N c 
corrections for the estimate of the d = 8 condensate contribu- 
tions; B(a s ) is the usual /^-function defined in Eq. (4) and A 2 is 
the square of the tachyonic gluon mass. 

• U c (q 2 )from MSW model 

The two-point correlator has been evaluated in the MSW model 
by [41, 42] in the Minkowski space and for a positive dilaton. 
Using the substitution in Eq. (23), we can transform the previ- 
ous result to the Euclidian space (q 2 < 0): 



IW) 



(d 2 ) 2 



■Ac 2 q 2 



log — 

V 1 



l[\2B(q 2 /c 2 g )-5] 



16 cl 



32 c£ 
'l5? 4 + 



(56) 



where one can also notice that the result is the same for a neg- 
ative dilaton. B(q 2 ) is a function of l/q 2 like in Eq. (25) but 
constant with respect to the variable z- In the MSW model, one 
takes B=0. In the HW model, one again finds that Yldq 2 ) does 
not have any power corrections like in the vector and scalar qq 
channels. 

• SVZ ffi \ I q 2 -expansion versus MSW model 

We compare the AdS/QCD expression in the Euclidian space 

with the one from the SVZ-expansion. 

- One obtains by matching the leading order PT expressions: 



8k 



B(a s ) 



where B(a s ) is defined in Eq. (4). 

- Equating the quadratic terms, one obtains: 

4cl = A 2 . 



(57) 



(58) 



- Equating the d > 4 terms, one finds in the MSW model that 
all of these terms have opposite signs with the ones in the SVZ- 
expansion. Such problems can be evaded by assuming that B + 
and by working within the expansion in Eq. (25) like done in 
[38]. However, the approach becomes obviously less predictive 
and would imply modifications in the Soft Wall recipe. 

6. Tachyonic gluon mass and hadron mass scale hierarchy 

• Tachyonic gluon and hadron spectra 

The parameter c? > which is the inverse of a squared length 
and which parametrizes the strength of a positive dilaton field 
cD,(z) = c 2 z 2 is assumed in the existing literature to be universal 
for a given model: 



2 2 2 
c — c — c 



(59) 



Therefore, from the mass relation for the n'''-radial excitation in 
the Minkowski space and [13, 38, 41]: 

fl - 4c 2 (n+l), 



M; 



Mi 



M, 



= 4< n+- , 



4c^(n + 2) 



(60) 



it is usual to deduce the ratios of the ground state masses (« 
0): 



3 

Ml * -M 2 n , M L r * 2M L n . 
2 



(61) 



One can also write the analogue of the previous results for the 
case of a negative dilaton field 0,(z) = -|c?|z 2 . One obtains 
[17]: 



M 



Pn 



4\c 2 ,\ (n + 1) 



M s„ = 4|c^n+-| 

M 2 Cn = 4\c 2 g \(n+l) 



1 



(62) 



If one assumes the equality in Eq. (59), one would obtain: 



M 2 



M* G *M* p 



(63) 



Ms is relatively low compared to the QCD spectral sum rules 
(QSSR) results for a qq meson mass Ms, ss 1 GeV [2], while 
the one of Mq is consistent with the one from a subtracted glu- 
onium sum rule: M Gl =a 1 GeV [40, 43, 44]. 
However, from our previous analysis in this paper and the re- 
lations between c, and the tachyonic gluon mass, one obtains 
instead of Eq. (59) the new relation for a negative dilaton: 



(64) 



indicating that e* is not necessary universal contrary to the 
tachyonic gluon mass A 2 . As shown in previous analysis, the 
value of the tachyonic gluon mass from phenomenological fits 
in Eq. (7) used in the holographic models leads to an un- 
derestimate of the QCD vacuum condensates. The value of 
A 2 which leads to the correct phenomenological values of the 
gluon (a s G 2 ) and four-quark (mm) 2 condensates ranges from the 
ones in Eqs. (38), (40) and (54): 



(*)• 



-(0.15 ~ 0.45) GeV 2 , 



(65) 



which is at the limit of the upper value obtained from phe- 
nomenological fit in Eq. (8). In the following discussion, we 
shall use the compromise value: 



(*)■ 



-(0.12 ~ 0.14) GeV 2 



(66) 



from the upper value in Eq. (8). From Eqs. (62) and (64), one 
can deduce in units of GeV: 



Mn 



0.5 , 



(67) 



which implies: 



M s ~ V2M p * 0.7 , M G * y/n/2a s M p * 0.8 , (68) 

where we have used (a s /7r) w — 1/ySi - 2/11 in the large N c - 
limit. Taking as a final estimate the range of values obtained in 
Eq. (61) from a positive and in Eq. (68) from a negative dilaton, 
one can deduce: 



Mi 

-4 ~ 1-5 

M 2 



2.0. 



M 2 
M 2 



2as s 



2.0 ~ 2.8 . 



(69) 



Assuming that the systematic uncertainties of the approach tend 
to cancel for the ratios of masses, we shall use the experimental 
value of M p = 0.778 GeV in order to deduce the final conser- 
vative results: 

M s ~ (0.95 ~ 1.1) GeV , M G * (1.1 ~ 1.3) GeV .(70) 

The masses of the vector and scalar mesons agree within 20% 
with the experimental value [p(0.78) and arj(0.98)] which is the 
expected accuracy of these leading large N c -limit holographic 
models. The predicted mass of the gluonium is in the range 
predicted by the subtracted (M G| 1 GeV) and unsubstracted 



(M G , =± 1.5 GeV) QCD gluonium spectral sum rules [40, 43, 
44] where the presence of two gluonia states are necessary for 
a consistency of the two sum rules \G\ is expected to couple 
strongly to pion pairs and might be identified to the <x(500) [40, 
43, 45] as indicated from recent fits of jui and yy scatterings 
data [46], while G2 is the "real glueball" which decays into the 
77' 77, ...U(l) channels]. The value of M G is consistent with the 
one from instanton sum rule [42] and with the lattice results of 
about 1 GeV including dynamical quarks [47] but lower than 
the quenched result of 1.6 GeV [48] 6 . 

• Hadronic mass scale 

Like in the QCD spectral sum rule analysis [3], one can also use 
c, for fixing the scale of the analyzed hadronic channel which 
is not necessary universal 7 . Following [3], one can deduce the 
hadronic scale from Eq. (64): 



2|c,r«4A^*4.2(a I /w)|^ 
4|cJ 2 * \A 2 \ * 5.5A 2 , 



(71) 



where we have used (a jn) « 2/11. The sizes of these scales 
are similar to the optimization scale appearing in the QCD spec- 
tral analysis of each channels but they are not necessary the 
hadron masses. Our results indicate the scale hierarchy: 

A p < A s < A G , (72) 

as expected in [3, 39]. 



7. Summary and conclusions 

We have systematically studied two-point functions built 
from light quark currents and gluon fields using the SVZ 
ffi l/q 2 expansion and compared it with the predictions from 
some QCD holographic models, namely the Minimal Soft-Wall 
Model (MSW model), the Gauge/String Dual Model (GSD 
model) and the Hard- Wall Model (HW model). These two dif- 
ferent approaches are expected to give a good parametrization 
of the non-perturbative aspect of QCD at moderate values of q 2 
in the Euclidian region. We found that in most of the channels 
studied, the HW model shows a large departure from the SVZ 
idea as it does not have any power corrections. We also notice 
the equivalence between the MSW model and the GSD model 
analyzed by Ref. [8-10] from a constraint among the length pa- 
rameters cj in the two models which is also obtained from their 
explicit relation with the tachyonic gluon mass. 
The quadratic 1 /g 2 -term introduced in [3-5] in order to quan- 
tify the UV renormalon contributions and large order of the 
PT series appeared naturally in all examples studied within the 
MSW model and the GSD model with a sign and size consistent 
with a tachyonic guon obtained phenomenologically [3, 15]. 



6 For a recent review on the glueball status, see e.g. [49]. Some other 
attempts to estimate the glueball masses in some holographic models can be 
found in [50]. 

7 These non-universality of c, can be also signaled by the difference of PT 
radiative corrections in each channel which has also been used in [5 1 ] for fixing 
the hadronic scale hierarchy. 
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Pursuing the matching of the two approaches for higher dimen- 
sion condensates, one notice that the correlator of the qq scalar 
current is described consistently term by term by the two ap- 
proaches until the d — 8 condensates contribution reached in 
the SVZ-expansion. In the vector channel, the MSW model 
fails to reproduce the SVZ coefficient of the d — 6 condensate 
contribution which vanishes in the MSW model due to some 
miraculous cancellation of its coefficients, while in the scalar 
gluonium channel, the MSW model cannot reproduce the sign 
of the Wilson coefficients of the d > 4 dimension condensates. 
In order to cure this difficulty, we have briefly discussed an ex- 
tension of the MSW model which is however less attractive than 
the MSW model due to its lesser predictivity. 
From the relation of the dilaton parameter c, to the tachyonic 
gluon mass, we have also deduced a constraint among these pa- 
rameters in different channels indicating that they are channel 
dependent but not universal for a given model. Then, we have 
deduced new values of the hadronic mass ratios where the pre- 
diction for the scalar gluonium mass is in the range given by the 
substracted and unsubtracted gluonium sum rules [40, 43, 44] 
and by lattice calculations with dynamical quarks [47]. We have 
also derived a hadronic hierarchy scale similar to the one ob- 
tained from QCD spectral sum rules analysis of different chan- 
nels [1-3, 39]. 
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